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ABSTRACT 

We study the three-dimensional and projected shapes of galaxy groups in the Sloan 
Digital Sky Survey Data Release 4, and examine the alignment between the orientation 
of the central galaxy and the spatial distribution of satellite galaxies. The projected 
ellipticity of a group is measured using the moments of the discrete distribution of its 
member galaxies. We infer the three-dimensional and projected axis ratios of their dark 
matter haloes by comparing the measured ellipticity distributions with those obtained 
from Monte Carlo simulations of projected, triaxial dark matter haloes with different 
axis ratios. We find that the halo shape has a strong dependence on the halo mass. 
While the haloes of low-mass groups are nearly spherical, those of massive groups tend 
to be prolate. For groups containing at least four members, the statistical distribution 
of their measured ellipticities does not have a strong dependence on the colors of their 
central galaxies. Our analysis further shows that the average three-dimensional axis 
ratio for haloes with 12 < log[M / {h~^ Mq)] < 15 is about 1 : 0.46 : 0.46, resulting in a 
projected axis ratio of ^ 0.77. Our results for the alignment between the orientation of 
the central galaxy of a group and the distribution of their satellite galaxies are in broad 
agreement with those obtained by Yang et al. The distribution of satellite galaxies 
preferentially aligns with the major axis of the central galaxy, with a clear dependence 
on both halo mass and galaxy colors. In particular, the alignment is stronger in more 
massive groups, and the strongest alignment is seen between red centrals and the 
distribution of red satellites. For groups with blue centrals, no significant alignment 
is detected. Finally, we examine how the observed alignment can be reproduced with 
the information about the halo axis ratios. The observed alignment signal can be 
reproduced if the angle between the major axis of the central galaxy and the projected 
major axis of the host halo has a Gaussian distribution with a mean of 0° and a 
dispersion of ^ 23°. This dispersion is larger for groups with blue centrals than those 
with red centrals. 

Key words: methods: statistical-galaxies: haloes- galaxies: structure-dark matter- 
large scale structure of universe. 



1 INTRODUCTION 

In the cold dark matter scenario, small dark matter haloes 
form first and grow subsequently to larger structures via ac- 
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cretion and merging processes. Such processes are generally 
anisotropic so that dark matter haloes are expected to be 
non- spherical. The orientations of dark matter haloes can 
be related to their surrounding structures, such as filaments 
and large-scale walls (e.g., Faltenbacher et al. 2002; Einasto 
et al. 2003; Avila-Reese et al. 2005; Hopkins, Bahcall & Bode 



2 Y. Wang et al. 



2005; Kasun & Evard 2005; Basilakos ct al. 2006; Altay et 
al. 2006; Aragon-Calvo et al. 2006; Maulbetsch et al. 2007; 
Ragone-Figueroa & Plionis 2007; Hahn et al. 2007a, 2007b). 
The last major merger of a dark matter halo may play an 
important role in determining its shape (van Haarlem & van 
de Weygaert 1993), although interactions between the gas 
and the dark matter components are also excepted to play 
a role (Kaaantzidis et aJ. 2004). Thus the shapes and orien- 
tations of dark matter haloes contain abundant information 
about their formation histories, which, in turn depend on 
the underlying cosmology (Lee 2006; Ho et al. 2006). On 
the other hand, the non-sphericity of dark matter haloes 
can also lead to systematic errors in cosmological studies 
(Sulkancn 1999; Wang & Fan 2006). Therefore it is of great 
importance to characterize the shapes of dark matter haloes, 
both observationally and theoretically. 

The shapes of clusters of galaxies can be probed using 
X-ray observations, studies of the Sunyaev-Zeldovich effect, 
and gravitational lensing (e.g., Zaroubi et al. 1998, 2001; 
Rehnsky 2000; Lee & Suto 2004; Wang & Fan 2004; De 
Filippis ct al. 2005; Scrcno ct al. 2006). However, one can 
also use the spatial distribution of satellite galaxies, since 
these are expected to be good tracers of the shapes of their 
host haloes. The large redshift surveys carried out in re- 
cent years, for example, the two-degree Field Galaxy redshift 
Survey (2dFGRS) and Sloan Digital Sky Survey (SDSS), 
have provided angular positions and redshifts for hundreds 
of thousands of galaxies, which allow detailed studies of the 
shapes of large samples and the shape dependence on rich- 
ness, multiplicity and dynamical evolution of groups and 
clusters. Early studies preferred the prolate shapes (Carter 
& Metcalfe 1980; Binggeli 1980; Phonis, Barrow & Frenk 
1991; Fasano et al. 1993; Cooray 2000; Basilakos, Plionis & 
Maddox 2000; Orlov, Petrova & Martynova 2001), but did 
not exclude the oblate solutions. Plionis et al.(2004) (here- 
after P04) estimated the shape distribution of UZC-SSRS2 
groups of galaxies by analyzing the spatial distribution of 
group members. They found that the prolate-like shape fits 
very well the cosmic structure on a large scale apart from 
the disc galaxy. Somewhat surprisingly, they also found that 
poor groups are more elongated than rich ones, results which 
are opposite to what is found in numerical simulations (eg., 
AUgood et al. 2006; Kasun & Evrard 2005). 

Somewhat surprisingly, these results are opposite to 
what is found with numerical simulations (Allgood et al. 
2006; Kasun & Evrard 2005). Recently, they estimated 
the average group morphology of the Percolation-Inferred 
Galaxy Group (2PIGG) and found that the prolate or, tri- 
axial with pronounced prolate shapes are the only acceptable 
morphological model (Plionis et al. 2006). 

In this paper, we use data from the SDSS to study the 
shapes of galaxy groups, under the assumption that the spa- 
tial distribution of group members traces the matter dis- 
tributions of their underlying dark matter haloes. Differing 
from previous studies, we run Monte Carlo simulations to 
generate the same number of member galaxies as in each 
of the observed groups under the assumption that haloes of 
the groups follow the triaxial model of Jing & Suto (2002) 
(hereafter JS02). We then make 2D projected distributions 
of these member galaxies, which serve as our Monte Carlo 
mock samples. We compare the ellipticity distribution ob- 



tained from the mock sample with the observed one to ex- 
tract information on the axis ratios of dark matter haloes. 

In the second part of this paper we revisit the align- 
ment between the spatial distribution of satellite galaxies in 
groups and the orientation of their central galaxies (hereafter 
'centrals'). It is important to assess accurately the align- 
ment of the dark matter halo because it holds important 
clues regarding the actual assembly history of dark matter 
haloes. Extensive studies with high resolution simulations 
have shown dark matter haloes have anisotropic distribu- 
tions of subhaloes that are aligned with their major axis 
(Knebe et al. 2004; Libeskind et al. 2005; Wang et al. 2005; 
Zentner et al. 2005). This anisotropy mainly owes to a pre- 
ferred direction of satellite accretion along large-scale fila- 
ments (Tormen 1997; Vitvitska et al. 2002; Aubert, Pichon 
& Colombi 2004; Knebe et al. 2004; Wang et al. 2005; Zent- 
ner et al. 2005). In addition, the tidal forces from the host 
halo may also induce new alignments (e.g., Ciotti & Dutta 
1994; Usami & Fujimoto 1997; Fleck & Kuhn 2003; Wang 
et al. 2007). Conversely, some nonlinear effects such as vi- 
olent relaxation and encounters can weaken the primordial 
alignment (e.g., Porciani et al. 2002). 

Comparing with the simulation studies, the observa- 
tional search for a possible alignment of central galaxies and 
their satellites has a long and confusing history. The first 
study of such an alignment was by Holmberg (1969), who 
found that satellites are preferentially located along the mi- 
nor axes of isolated disc galaxies. Holmberg's study was re- 
stricted to projected satellite-central distances of rj, < 50kpc. 
Subsequent studies, however, were unable to confirm this 
so-called 'Holmberg effect' (Hawley & Peebles 1975; Sharp, 
Lin & White 1979; MacGillivray et al. 1982). Zaritsky et 
al. (1997) studied the distribution of satellites around spiral 
hosts and were also unable to detect any significant align- 
ment for Vp < 200 kpc, but they found a preferred minor-axis 
alignment for 300 kpc <rp< 500 kpc. Our Milk Way and 
M31 have satellites that lie in great planes that are highly 
inclined to their discs (Lynden-Bell 1976, 1982; Majewski 
1994; Hartwick 1996, 2000; Kroupa et al. 2005; Koch & 
Grebel 2006; McConnachie & Irwin 2006; Metz et al. 2007). 
With large redshift surveys, such as 2dFGRS and SDSS, 
much larger samples of galgixy groups can be used to dis- 
cuss the alignment problem. Sales & Lambas (2004) used a 
set of 1498 host galaxies with 3079 satellites from the 2dF- 
GRS and found a large-scale ( alignment of the satellites 
along the host minor axes for 300 kpc < r-p < 500 kpc. Brain- 
erd (2005) studied a sample of isolated SDSS galaxies and 
found that the distribution of satellite galaxies is strongly 
aligned with the major axis of the disc host galaxy. Yang 
et al. (2006, hereafter Y06), using a galaxy group catalogue 
similar to the one used here, but based on the SDSS Data 
Release 2 (DR2), studied the alignment signal £is function of 
the colors of the central and satellite galaxies. They found 
that the alignment strength is strongest between red centrals 
and red satellites, while the satellite distribution in systems 
with a blue central galaxy is consistent with being isotropic. 
Y06 also found that the alignment strength is stronger in 
more massive haloes and at smaller projected radii from 
the central galaxy. These results have subsequently been 
confirmed by several independent studies (Donoso, O'Mill 
& Lambas 2006; Azzaro et al. 2007; Agustsson & Brain- 
erd 2006a, 2007). Using the same group catalogue as that 
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used here, Faltenbacher et al. (2007a) examined several ad- 
ditional alignment signals. They found that the orientations 
of red satellites are preferentially aligned radially in the di- 
rection of the brightest group galaxies (BGG). In addition, 
they found a weak but significant indication that the orien- 
tations of satellite gala^xies are directly aligned with that of 
their BGG. Gomparing with the earlier studies of the align- 
ment between brightest cluster galaxies (BCGs) and their 
parent clusters (Carter & Metcalfe 1980; Binggeli 1982; Stru- 
blc 1990; West 1994; Kirn ct al. 2001), Faltenbacher ct al. 
(2007a) have given more detailed results for the large sam- 
ples. These various detections of alignment between centrals 
and satellites have triggered a number of investigations into 
the connection between the shapes and orientations of dark 
matter haloes and their galaxy population, with the goal to 
improve our understanding of the formation of dark matter 
haloes and galax:ies (e.g. Agustsson & Brainerd 2006b, here- 
after, AB06; Kang et al. 2007, hereafter K07; Faltenbacher 
et al. 2007b; Brunino et al. 2007; Sales et al. 2007; Pereira 
et al. 2007). 

In addition to measuring the alignment signals from the 
SDSS observations for groups of different masses, we try to 
infer the correlation between the orientations of the central 
galaxy and that of its host halo, statistically but both from 
the observations. If the projected orientation of a central 
galaxy is perfectly aligned with the projected orientation of 
its host dark matter halo, assuming that satellite galaxies 
trace the matter distribution, the alignment signal between 
the distribution of satellite galaxies and the orientation of 
their central galaxy is strongest. By comparing the align- 
ment signals measured from the observations and the Monte 
Carlo samples, we can estimate the deviation (misalignment 
angle) of the orientation of central galaxy from the orienta- 
tion of its host dark matter halo. We use a Gaussian distri- 
bution function to quantify this misalignment angle. 

This paper is organized as follows. In §2, wc briefly 
describe the observational data used for this study. Sec- 
tion 3 presents our measurements of the intrinsic shapes of 
dark matter haloes, where the three-dimensional and two- 
dimensional axis ratios are determined using Monte Carlo 
simulations. Section 4 shows the alignment signal we mea- 
sured from the SDSS and its implication for the shape cor- 
relation between the central galaxies and the dark mat- 
ter haloes. Finally, §5 presents a summary and discussion. 
Throughout this paper we refer to the inferred shape from 
the satellite galaxy distribution as the shape of the group 
and the corresponding dark matter halo, and use the major- 
axis direction of the satellite distribution to indicate the 
orientation of the group and the corresponding dark matter 
halo. 



in the range 0.01 < z < 0.20 and with a redshift complete- 
ness C > 0.7. This sample of galaxies is used to construct 
three group samples: sample I, which only uses the 362356 
galaxies with measured redshifts from the SDSS, sample 
II which also includes 7091 galaxies with SDSS photome- 
try but with redshifts taken from alternative surveys, and 
sample III which includes an additional 38672 galaxies that 
lack a redshift due to fiber-collisions, but which we assign 
the redshift of its nearest neighbor (cf. Zehavi et al. 2002). 
The present analysis is based on sample II which consists of 
369447 galajcies distributed over 301237 groups with a sky 
coverage of 4514 deg^. Details of the group finder and the 
general properties of the groups can be found in Yang et al. 
(2007). 

In this paper, the central galaxy is defined to be the 
brightest galaxy in the group and other galaxies are satel- 
lites. We also take the most massive (in terms of stellar 
mass) group member as the central galaxy. As we have 
tested, the difference between these two definitions is too 
small to be noticed. The group masses are estimated using 
the ranking of group's characteristic luminosity, Lig.s, de- 
fined as the combined luminosity of all group members with 
"■^A/,. — 51og/i < —19.5. More details of the mass estima- 
tions can be found in Yang et al. (2007). Note that, the 
selected galaxy groups contain a small fraction of interlop- 
ers, i.e. false members assigned to a group. If the distribution 
of these interlopers is uncorrelated (or anti-correlated) with 
that of the true members of the group, our results on both 
the ellipticity and the alignment can be biased. According 
to Yang et al. (2007; 2005a) the average fraction of inter- 
lopers in the group is less than 20%. We have tested the 
effect of such fraction by assuming that the distribution of 
the interlopers is uncorrelated with the shape of of the group 
and is spherical, we find that the presence of the interlopers 
can decrease the ellipticity of the groups and the alignment 
signals by ^ 10%. 

Note that in these group catalogues survey edge effects 
have been taken into account (Yang et al. 2007). Only groups 
with /edge > 0.6 are selected, where 1 — /edge is the fraction 
of galaxies in a group that are missed due to the edge ef- 
fects. In order to obtain the ellipticity distribution of galaxy 
groups, we only use groups with at least four members (one 
central galaxy and at least three satellites), which results in 
a catalogue of 5184 groups. However, in studying the align- 
ment between satellite galaxies and the orientation of their 
centrals, wc enlarge our sample by using all groups with at 
least two members (one central and one satellite). This sam- 
ple gives a total of 62212 unique central-satellite pairs, many 
more than in the DR2 sample used by Y06. 



2 DATA 

The analysis presented in this paper is based on the SDSS 
DR4 galaxy group catalogue of Yang et al. (2007). This 
group catalogue is constructed applying the halo-based 
group finder of Yang et al. (2005a) to the New York Univer- 
sity Value-Added Galaxy Catalogue (NYU-VAGC; see Blan- 
ton et al. 2005), which is based on SDSS DR4 (Adelman- 
McCarthy et al. 2006). Prom this catalogue Yang et al. se- 
lected all galaxies in the Main Galaxy Sample with redshifts 



3 THE INTRINSIC SHAPE OF THE DARK 
MATTER HALO 

3.1 Methodology 

We now describe how we use the satellite distribution to 
determine the ellipticity distribution of galaxy groups and 
their corresponding dark matter haloes. The observed satel- 
hte distribution in a group suffers from severe discreteness 
effects. In particular, since each group contains only a small 
number of galaxies, there is a high level of Poisson noise in 
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the determination of the eUipticity based on its galaxy dis- 
tribution. Thus the eUipticity directly measured can only be 
used as a rough indicator of the underlying, true eUiptic- 
ity. We will use mock samples to quantify how the observed 
eUipticity distribution is related to the real distribution. As- 
suming that the distribution of satellite galaxies in a group 
traces the mass distribution in the corresponding dark mat- 
ter halo, we can infer, in a statistical sense, the shapes of 
dark matter haloes from the observed distribution of the 
group eUipticities. In order to obtain the principal axes and 
the orientation of a group projected on the sky, we define 
the inertia tensor as 



(1) 



where {xi^„,Xj^„) are the projected coordinates (with the 
central galaxy at the origin) of the n"* satellite galaxy. The 
semi-major and semi-minor axes of the ellipse. La and Lt 
(two roots of the following equation) , can be derived by solv- 
ing the equation 

r 2 

= 0. (2) 
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The direction of the major axis is given by the eigenvector 
r — [1, (L^ — X\i) / Xi2\, while the eUipticity, e, and the axis 
ratio, r], are 



e = l-Lt/La, and r] = Li/La- 



(3) 



Throughout this paper, we use the eUipticity e and axis ratio 
r] to refer to the quantities measured directly from the data. 
The inferred shapes of the underlying dark matter haloes 
are described either by their 2D or 3D axis ratios. 

As mentioned above, to quantify the true shapes of 
the dark matter haloes associated with the galaxy groups, 
one needs to construct mock samples that include the same 
selection effects. We construct Monte Carlo SDSS DR4 
group catalogues as follows. First, we determine the num- 
ber of satellites for each of the SDSS groups. Second, we 
re-distribute these satellites according to a spherical NFW 
profile (Navarro et al. 1996, 1997) or a triaxial profile by 
JS02. Third, we project the three dimensional distribution of 
satellite galaxies onto a two dimensional plane and measure 
the eUipticity, emc, for each Monte-Carlo group. Such Monte 
Carlo approach was firstly introduced by Basilakos, Plionis 
& Maddox (2000) to recover the projected cluster eUipticity 
distribution and the true projected eUipticity distribution 
taking into account the background and discreetness efi'ects 
(see also Plionis et al. 2006). Thirty Monte-Carlo realiza- 
tions are generated, and we estimate the average and scat- 
ter of the eUipticity distribution using these thirty Monte 
Carlo samples. Some small groups with masses smaller than 
lO"-® /i"^M0, where SDSS DR4 group catalogue does not 
provide mass estimates, are removed from our sample. 

The first model for the mass profile used in our Monte 
Carlo simulations is a spherical NFW density profile, 

po5c 



PNPw(r) = 



(r/r,)(l + r/r,)2 



(4) 



where po is the average mass density of the universe, is 
a scale radius and 5^ = 200cV3[ln(l + c) - c/{l + c)]. The 
concentration parameter c is defined as c = r2oo/''3, with 
'"200 the radius within which the mean density is 200 times 



the average mass density of the universe. This model is fully 
determined for a halo of a given mass (or, equivalently, r2oo) 
once the concentration parameter is given. The concentra- 
tion parameter depends on the halo mass M and redshift z, 
for which we use the model of Bullock et al. (2001): 



c{M,z) 
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M 
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where c* = 8, as is appropriate for the ACDM model. 

The other density profile we use is that proposed by 
JS02. Using high-resolution numerical simulations, JS02 
proposed an NFW-like triaxial density profile for dark mat- 
ter haloes, which has the form 



P(^).7 



poSc 



{R/Ro)"'{l + R/Roy 



(6) 



where R = a(a;^/a^-f y^/fo^+z^/c^)^/^, and a > b > c are the 
lengths of the three principal semi-axes. For the value of q, it 
is found that both a = 1 and 1.5 can provide a good fit to the 
simulated profiles. Detailed comparisons showed that a = 1 
is slightly better for cluster-scale haloes, while a — 1.5 gives 
better fit for galactic haloes (JS02). Therefore, we adopt 
Q — 1.5 for groups with masses smaller than 10^"^ h~^MQ, 
and Of = 1.0 for more massive groups (In fact, the results 
are not sensitive to the value of a adopted. We have also 
adopted a = 1.0 for all groups, their differences are tiny.). 



3.2 EUipticity distribution 

Figs. [1] and [2] show the probability distributions of the eUip- 
ticity obtained from the SDSS DR4 (dotted histogram) and 
from the Monte Carlo simulations (solid histogram with er- 
ror bars). In Fig.[T] we assume that the distribution of satel- 
lite galaxies follows a spherical NFW profile, while in Fig. [2] 
we assume that the distribution of satellite galaxies follows 
the JS02 model. For the JS02 triaxial distribution, we have 
adopted the model parameters (axis ratios and concentra- 
tions) given in JS02. In both cases we randomly select a 
line-of-sight direction, project the three-dimensional satel- 
lite distribution onto a two-dimensional plane, and measure 
the eUipticity distribution using the method outlined above. 
In each plot, the parameter represents the lower limit on 
the number of galaxies in each group (including the central 
galaxy). Note that the eUipticity distribution is strongly de- 
pendent on this lower limit, with poorer groups being more 
elongated. The mean eUipticity increases from ~ 0.40 for 
groups with richness A > 10 to ~ 0.54 for groups with 
A' > 4. However, this does not mean that the true halo shape 
is more elongated for poorer groups; the trend is largely a 
result of the discrete sampling. For example, in the extreme 
case where only one satellite galaxy is observed, the mea- 
sured eUipticity will always be unity regardless of the shape 
of the underlying dark matter halo. Thus, the shape of the 
underlying halo shape can only be probed in an indirect way, 
i.e., by comparing the observed distribution with that of the 
Monte-Carlo samples. 

From Fig. [T] and [51 it appears that the triaxial model 
of JS02 fits the data better relative to the spherical NFW 
model. 
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Figure 1. Probability distributions of the ellipticity in the SDSS 
DR4 (dotted histogram) and Monte Carlo re-samphngs (sohd his- 
togram with error bars). Here we assume that the distribution of 
satellite galaxies in groups follows a spherical NFW profile. 
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Figure 2. Same as Fig. [T] but here we assume that the satel- 
lite galaxies are distributed in groups according to JS02 triaxial 
model. 



3.3 The three-dimensional and projected shapes 
of dark matter haloes 

Numerous studies have used numerical A''-body simulations 
to probe the non-spherical shapes of dark matter haloes as 
traced by dark matter particles (e.g., JS02; Kazantzidis et 
al. 2004), or by subhaloes (e.g., Diemand, Moore & Stadel 
2004). The subhaloes are more closely associated with satel- 
lite galaxies, which have been found to be biased tracers of 
the mass distribution. There is a negative spatial bias at the 
center, and a corresponding positive velocity bias. It is still 
unclear whether this reflects numerical artefact (i.e., over- 
merging), or whether this is real. Yang et al. (2005b) studied 
the satellite distributions in the 2dFGRS groups and found 
evidence that the number density distribution of satellites 
is less concentrated than expected dark matter. Here we 
probe the mean values of axis ratios of dark matter haloes 
as traced by SDSS galaxies. For this purpose, we first divide 
the SDSS DR4 groups into subsamples according to their 
halo masses, and measure the corresponding ellipticity dis- 



Table 1. Best-fit parameters for groups in different halo mass 
bins. From top to bottom, values are listed for groups with dif- 
ferent members A'^ > 4, 6, 8, and 10, respectively. 



Mass bin 3D axis ratios 2D axis ratios 
(I) (I) (II) (11) 

12 < M' < 13° 1:0.96:0.85 22.58 0.96^0.02 46.16 

13 < M' < 14 1:0.80:0.72 29.93 O.SStoM ^^'^O 

14 < M' < 15 1:0.88:0.44 11.86 0.69t° ol 19.13 
12 < M' < 15 1:0.46:0.46 67.63 0.77to.o? 56.02 



12 < M' < 13 1:0.96:0.95 12.47 O.86I0 ^4 20.49 

13 < M' < 14 1:0.96:0.63 10.75 0.76lg |j} 16.29 

14 < M' < 15 1:0.35:0.35 17.78 0.68to.ol 52.10 
12 < M' < 15 1:0.51:0.48 13.17 0.72to o} 14.63 



12 < M' < 13* 

13 < M' < 14 1:0.87:0.53 5.36 0.72lg gi 10.40 

14 < M' < 15 1:0.33:0.33 33.19 O.lOtoxa 89-14 
12 < M' < 15 1:0.43:0.43 17.05 0.72t^;[5} 29.53 



12 < M' < 13 

13 < M' < 14 1:0.84:0.50 7.68 0.72lg gi 12.63 

14 < M' < 15 1:0.61:0.36 34.84 0.68to °i 60.24 
12 < M' < 15 1:0.39:0.39 23.04 0.7lt° oi 48.66 



" M' = log[M/(/i-iM0)]. 

^ The axis ratios in the 12 < M' < 13 mass bin with N > 8 and 
Af > 10 are absent because our group sample contains too few 
of these rich groups for a reliable measurement of the ellipticity 
distribution. 

tribution as traced by the satellite galaxies. Then, using the 
model of JS02 with given axis ratios, we construct thirty 
realizations of Monte Carlo simulations and measure the 
corresponding ellipticity distributions. By changing model 
parameters (i.e. axis ratios) so that the predicted elliptic- 
ity distributions match the observed one, we determine the 
underlying axis ratios of dark matter haloes. Here the com- 
parison between model predictions and observation is done 
in terms of a x^, defined as 

2 _ {{fi(MC)) — .fi(obs))^ , - 

where Ni, = 10 denotes the bin number of the ellipticity 
distribution. {fi(MC)) a-nd o-'^{fi(MC)) a-re, respectively, the 
average amplitude and 1-a deviation of the ellipticity dis- 
tributions obtained from the thirty realizations of Monte 
Carlo simulations, while fi(obs) is the amplitude of the el- 
lipticity distribution obtained from the SDSS groups. Note 
that for each set of axis ratios, {fi(MC) ) changes slightly. 
However, as we have tested, using constant cr^ifi^MC)) does 
not have a significant impact on our measurement of the 
best-fit axis ratios. Thus we use o-^ {fi(MC)) computed from 
the thirty Monte Carlo simulations and estimate the best-fit 
axial ratios by minimizing the x^- Here the minimization is 
performed on the regular 100 x 100 (3D) or 100 (2D) axis 
ratio grids. 

In Fig. [3] we compare the ellipticity distribution for 
groups in different mass bins (dotted histograms) with 
those of the best-fit Monte-Carlo simulation (solid his- 
tograms with errorbars). The corresponding best-fit, three- 
dimensional (3D) axis ratios for the dark matter haloes are 
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12<log[M/h"'MjS13 1 3<log[M/h"'MjS 1 4 1 4<log[M/h"'M JS 1 5 1 2<log[M/h"'MjS1 5 




0.0 0.4 0.8 0.0 0.4 0.8 0.0 0.4 0.8 0.0 0.4 0.8 

ellipticity ellipticity ellipticity ellipticity 

Figure 3. Similar to Fig. [T] but for groups in different mass bins as indicated on top of the panels. In this plot, the Monte Carlo 
simulations are performed according to the JS02 profile with the best-fit three-dimensional axis ratios. The line-of-sight (projection) 
direction is random. 



listed in the second column of Table [T] As one can see, the 
3D axis ratios recovered are different for groups with differ- 
ent richness. The results indicate that less massive haloes are 
more spherical, while more massive haloes tend to be more 
prolate. These results are in good agreement with Paz et al. 
(2006; hereafter P06), who studied the shapes of dark mat- 
ter haloes, both projected and in 3D, using simulations and 
groups constructed from the 2-degree Field Galaxy Redshift 
Survey (2dFGRS; CoUess et al. 2001) and the SDSS DR3. 

In addition to the three-dimensional axis ratios, we can 
also look into the projected two-dimensional (2D) axis ra- 
tios. Assuming a projected major-to-minor axis ratio, we re- 
sample the sky positions of satellite galaxies in each group in 
a Monte Carlo way, and measure the corresponding elliptic- 
ity distribution for groups of different masses and richness. 
Here again, by changing the 2D axis ratio, we can obtain the 
one that best matches the observed ellipticity distribution of 
groups. The results are shown in Fig.|4l and the correspond- 
ing best-fit, 2D axis ratios are listed in the fourth column of 
Table [1] Once again less massive groups are found be more 
spherical. Note that groups in the same mass bin, but with 
different richness, have only slightly different 2D axis ratios. 

The values presented in Table [T]for the best fit mod- 
els are quite large, especially for groups with halo masses 



10^^ - 10^'' /i-^Mq in the 3D cases, and in a few 2D cases. 
Possible reasons for such large are (i) the assumption of 
a constant axis ratio is not accurate to describe the shape 
of the groups, and (ii) the shapes of observed groups do not 
follow a Gaussian distribution. Unfortunately, without bet- 
ter knowledge about these issues, it is difficult to come up 
with a better model. 

Observationally, only gravitational lensing data can di- 
rectly probe the 2D, projected mass distribution of dark 
matter haloes. However, at the present time, such observa- 
tions are limited in both quality and quantity (see Hoek- 
stra et al. 2004; Mandelbaum et al. 2006). For clusters of 
galaxies, additional data is available from X-ray data and 
from studies of the Sunyaev-Zeldovich-effect. These data 
yield constraints on the shapes of cluster haloes from their 
gas distribution. Mohr et al. (1995) used a sample of 57 
X-ray clusters observed by the Einstein telescope and ob- 
tained a mean 2D axis ratio, {r7)=0.80, and a dispersion 
Gr, = 0.12. With the JSQ2 model, Wang & Fan (2004) studied 
the Sunyaev-Zeldovich effect and X-ray surface brightness 
profiles for clusters of galaxies. Based on a sample of clusters 
with masses above Mum = lO^^fe~^M0, they found that the 
average axis ratios is {rf) ~ 0.84. Sereno et al. (2006) used 
a sample of 25 X-ray selected clusters observed by Chan- 
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Figure 4. Similar to Fig. |3] but for the best-fit two-dimensional axis ratios. 



dra and XMM-Newton, and obtained a mean projected axis 
ratio (ri) = 0.80 ± 0.02. Recently, Flores et al. (2007) pre- 
dicted that the mean axis ratio is (rj) — 0.82 and a scatter 
an — 0.09 using a simple analytical model. Note, though, 
that all these measurements are for the axis ratios of the hot 
intra-cluster gas. Lee & Suto (2003) presented an analytical 
expression that links the ellipticity of the gas to that of the 
dark matter halo, assuming that the intra-cluster gas is in 
hydrostatic equilibrium. Although this relation is obtained 
in 3D, we assume that it is also valid for the 2D distributions. 
By converting the axis ratio of the gas distribution into the 
dark matter distribution using the results in Fig. 3 of Lee & 
Suto (2003), we infer 2D axis ratios for dark matter haloes 
of 0.6±0.1, in good agreement with our measurements (0.69 
for groups with halo masses 10^* — 10^® h~^MQ). 



Yang et al. (2008) we separate galaxies into red and blue 
population using the following dividing curve, 

"■^(9 - r)o = 1.022 - 0.0651 * Mr,23 - 0.00311 * M^^s , (8) 

where M,.,23 = '''^M,. - 5logh + 23, and "'^M^ - 5 log /i is 
the absolute magnitude K + E corrected to z — 0.1 using 
the method described in Blanton et al. (2003). We define 
galaxies with '^'^{g — r) >"'^ {g — r)o as red galaxies and the 
rest as blue galaxies. Here '^'^{g — r) is the color in the SDSS 
g and r bands K + E corrected to z = 0.1. In Fig. [5] we 
compare the ellipticity distribution for groups with red and 
blue central galaxies. The results indicate that there is no 
significant difference between the shape of groups with blue 
or red central galaxies. 



3.4 Dependence of halo shape on the color of 
central galaxies 

The color-magnitude relation of galaxies is found to have 
a bi-modal distribution, consisting of a red 'sequence' and 
a blue 'cloud' (e.g. Baldry et al.2004; Li et al. 2006). In 
this subsection we test whether haloes that host red or blue 
central galaxies have different shapes (the motivation behind 
this will become clear in Following Li et al. (2006) and 



4 ALIGNMENT BETWEEN CENTRAL 

GALAXIES AND DARK MATTER HALOES 

4.1 Quantifying the alignment 

In order to quantify the distribution of satellite galaxies in 
groups relative to the orientations of their central galax- 
ies we follow Brainerd (2005) and compute the distribution 
function, P{0), where 9 is the angle between the major axis 
of the central group galaxy and the direction of a satellite 
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Figure 5. Probability distribution of the ellipticity for groups 
with different central galaxies: red (solid line) v.s. blue (dotted 
line) . 



relative to the central galaxy. The angle 6 is constrained in 
the range 0° < 6* < 90°, where = 0°(90°) implies that 
the satellite lies along the major (minor) axis of the cen- 
tral galaxy. The orientation of the central galaxy is based 
on the isophotal position angle in the r band, as given in 
the SDSS-DR4 (Adelman- McCarthy et al. 2006). We have 
checked the distribution of these position angles and found 
it to be isotropic. 

For a given set of central and satellite galaxies, we first 
count the total number of central-satellite pairs, N{6), for a 
number of bins in 9. Next, we construct 100 random samples 
in which we randomize the orientations of all central galax- 
ies, and compute {Nii{9)}, the average number of central- 
satellite pairs as function of 6. Note that this ensures that 
the random samples have exactly the same selection effects 
as the real sample, so that any significant difference between 
N{9) and N[i{9) reflects a genuine alignment between the 
orientations of the central galaxies and the distributions of 
their corresponding satellite galaxies. 

To quantify the strength of any possible alignment we 
follow Y06 and define the distribution of normalized pair 
counts: 



/pairs {9) — 



N{9) 
{Nr{9)}- 



(9) 



Note that in the absence of any alignment, /pairs(S) ~ 1, 
while /pairs (S) > 1 at small 6 implies a satellite distribution 
with a preferred alignment along the major axis of their cen- 
tral galaxy. We use an{9) / {N-r{6)), where ctr is the standard 
deviation of N-r{9) obtained from the 100 random samples, 
to assess the significance of the deviation of /pairs (S) from 
unity. In addition to this normalized pair count, we also com- 
pute the average angle {6} . In the absence of any alignment 
(6) = 45°, however, {9} = 45° does not mean an isotropic 
distribution. Major and minor axis alignments are charac- 
terized by (6) < 45° and (9) > 45°, respectively. The sig- 
nificance of any alignment can be expressed in terms of ag, 
the variance in {6}r, which is obtained from the 100 random 
samples. 

Fig.[6]shows /pairs (^) for all central-satellite pairs (solid 
line) in our SDSS group catalogue. Clearly /pairs (^) > 1 
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Figure 6. The normalized probability distribution, /pairs(^)i of 
the angle between the major axis of the central galaxy and the 
direction connecting the satellite galaxy and the central galaxy. 
Results are measured for groups with at least two members (solid 
line) and four members (dotted line). 



for small 9, indicating that satellite galaxies are distributed 
preferentially along the major axis of their central galaxy. 
This is also evident from the fact that {9} = 42.46° ± 0.12°, 
which deviates from the case of no alignment (i.e., (9) — 45°) 
by almost 21a\ For comparison, using the a group cat- 
alogue constructed from the SDSS DR2 data by Wein- 
mann et al. (2006), similar to that used here, Y06 found 
(9) = 42.2° ± 0.2°, in excellent agreement with the results 
presented here. Note, though, that the statistical error pre- 
sented here is much smaller, due to the larger group cata- 
logue used. The existence of alignment owns in part to the 
non-spherical distribution of the satellite galaxies in dark 
matter haloes (e.g., Zentner et al. 2005; Kang et al. 2005, 
2007; Libeskind et al. 2005; AB06), which has been used in 
the previous section to probe the overall shapes of the dark 
matter mass distribution. For comparison, the dotted lines 
in Fig. [6] show /pairs (^) for groups with at least four mem- 
bers. The resulting alignment signal is slightly stronger than 
for the full sample (which includes all central-satellite pairs 
in groups with at least two members). This is consistent 
with the fact that (i) groups with more satellites are more 
massive and (ii) more massive groups are less spherical (see 
also gT2] below). 



4.1.1 Dependence on galaxy color 

In order to study how the alignment depends on various 
properties of the central and satellite galaxies, we follow Y06 
and divide our sample of central-satellite pairs into different 
subsamples. The upper panels of Fig. [7] show the alignment 
signals obtained for blue and red satellites, while the lower 
panels show the results for blue and red centrals. As one 
can see, there is a strong dependence on the colors of both 
the centrals and satellites. In particular, groups with red 
centrals and red satellites show a stronger alignment than 
those with blue centrals and blue satellites, in good agree- 
ment with previous studies (Y06; Azzaro et al. 2007; Agusts- 
son & Brainerd 2007). As pointed out by K07, groups with 
blue centrals tend to have slightly more interlopers, which 
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Figure 7. Same as Fig. |6] but for different subsamples of central- 
satellite pairs, separated according to the "'^(g — r) colors either 
for satellite galaxies (upper two panels) and for cental galaxies 
(lower two panels). 
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Figure 9. Similar to Fig.|6] but for central-satellite pairs in haloes 
of different masses as indicated on top of the panels. Results are 
shown separately for all, red and blue central galaxies in the up- 
per, middle and lower panels. 
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Figure 8. Same as Fig. except that here we split the sample 
according to different combinations of the colors for both central 
and the satellite galaxies, as indicated. 

may cause an arteficial reduction of their measured align- 
ment. However, even after this is corrected for, groups with 
blue centrals still show a weak alignment. 

In Fig. [8] we show the alignment , /pairs (6*), for the four 
color combinations of centrals and satellites. As one can 
see, pairs between blue centrals and blue satellites do not 
show any alignment, while pairs between red centrals and 
red satellites show the strongest alignment. Pairs between 
red centrals and blue satellites and pairs between blue cen- 
trals and red satellites show alignment with intermediate 
strength. All these findings are in excellent agreement with, 
and more significant than, those obtained in Y06. 

4.1.2 Dependence on halo mass 

Fig. [9] shows the alignment measure for groups of different 
halo masses. From the upper panels one can see that the 
alignment is stronger for more massive groups. We also ex- 
amine the mass dependence of the alignment separately for 



groups with blue and red centrals, the results of which are 
shown in the middle and lower panels of Fig. [9j respectively. 
The alignment is quite different for blue and red centrals. 
For all halo masses probed here, red centrals reveal a much 
stronger alignment with their satellite galaxies than blue 
centrals. In fact, except for the most massive haloes with 
14 < logAf/ h~^MQ < 15, we find no significant alignment 
signal of satellites with blue centrals. These results are again 
in good agreement with the findings of Y06 based on galaxy 
groups in the SDSS DR2 (see also Agustsson & Brainerd 
2007). 

4.2 The alignment of central galaxies with their 
host haloes 

As shown above, central galaxies are aligned with the distri- 
bution of their satellite galaxies, and different systems show 
different levels of alignment. To produce such an alignment 
signal requires that (i) dark matter haloes are elongated, 
(ii) satellite galaxies trace the dark matter, at least to some 
extent, and (iii) the orientation of the central galaxy is corre- 
lated with the orientation of its dark matter halo. Now that 
we have obtained constraints on both the shapes of dark 
matter haloes, we can use the observed alignment signals to 
constrain the correlation between the orientations of central 
galaxies and their dark matter haloes. 

Previously, such correlations have been constrained us- 
ing galaxy catalogues constructed from semi-analytical mod- 
els (SAM) for galaxy formation (see AB06 and K07). In K07, 
the satellite distributions are modelled using the locations of 
subhaloes (which are thought to host the satellite galaxies) 
inside larger dark matter haloes in a numerical simulation. 
In general, subhaloes are found to be accurate tracers of the 
shapes of the dark matter distributions of their host haloes. 
However, since the SAM do not predict the orientations of 
the central galaxies, a number of simple assumptions have 
been made so far, where the minor axis of the central is per- 
fectly aligned with (i) the major axis of the inertia tensor of 
the host halo, (ii) the minor axis of the inertia tensor of the 
host halo, (iii) the intermediate axis of the inertia tensor of 
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Figure 10. The alignment signal /pairs C^') obtained by consider 
the misalignment-angle between the major axis of the dark mat- 
ter halo and the major axis of the central group galaxy. Upper 
panels: the alignment signals by assuming that the direction of 
central galaxy is same as the direction of the projected major axis 
of the dark matter halo (dashed line), comparing to the observed 
alignment signal (solid line). Lower panels: The comparisons be- 
tween the observed alignment signal (solid line) and the best fit- 
ting alignment signal by considering a misalignment between the 
projected major axis of the dark matter halo and the major axis 
of the central galaxy (dashed line). 

the host halo, or (iv) the angular momentum vector of the 
host halo. Both AB06 (using the alignment signal measured 
from isolated host-satellite systems) and K07 (using central- 
satellite pairs in galaxy groups) found that model (i) predicts 
an alignment signal that is much too strong compared with 
the data, model (ii) predicts a strong Holmberg effect, con- 
trary to what is seen, and model (iii) predicts almost no 
alignment. Only model (iv) results in alignment signals that 
are in agreement with observations. Although this does not 
give definite proof that the minor axis of the central galaxy is 
perfectly aligned with the angular momentum vector of the 
halo, it does provide a possible understanding for the origin 
of a correlation between the orientation a central galaxy and 
that of its host halo. 

However, it is important to keep in mind that the SAM 
not necessarily predict the correct spatial distributions of 
satellite galaxies. Therefore, in this paper we focus on what 
can be inferred from the data alone. In particular, we try to 
mfer, based purely on the data presented here, to what ex- 
tent the projected orientations of central galaxies are aligned 
with those of their dark matter haloes. Similar to what 
we have done to extract the three-dimensional and two- 
dimensional axis ratios in Section |3] we compare the data 
to Monte Carlo simulations. We start from the groups with 
at least four members, and assume that the projected orien- 
tations of the central galaxies are perfectly aligned with the 
orientation of the projected satellite distribution. Next, us- 
ing a Monte Carlo method, we distribute the satellite galax- 
ies according to the projected two-dimensional axis ratios 
listed in the fourth column of Table [1] for the group sample 
with A'' > 4. Finally, we measure the various alignment sig- 
nals for this Monte Carlo simulation using the same method 
as described above. The upper panels of Fig. [10] show the 
alignment signals thus obtained (open symbols) for groups 



of difi'erent masses, as indicated at the top of each panel. 
For comparison, we also plot the data as asterisks. Clearly, 
this model predicts an alignment signal that is stronger than 
observed, especially for the more massive groups. 

We can suppress the strength of the alignment signals in 
the model by taking scatter (i.e., 'random' deviations from 
perfect alignment) into account. To that extent we assume 
that the misalignment angle between the projected orienta- 
tion of the central galaxy and the major axis of its host halo 
can be described by a Gaussian distribution: 

P(...)^^exp(-||), (10) 

with <j the standard deviation of the distribution. Note that 
this distribution for the misalignment angle Omia is symmet- 
ric and centered around zero. 

Using such a Gaussian distribution, we fit the various 
alignment signals obtained from our SDSS group catalogue, 
treating cr as a free parameter in each separate case. Here 
again we use the minimum fit, similar to Eq.[7]but for the 
angular distribution of central-satellite pairs (Eq. The 
parameter o is constrained in the range 0° — 90° , and on the 
regular 900 grids. For each o, the misalignment angles Smis 
are generated according to the Gaussian distributions. We 
measure the alignment signals in Monte Carlo mock samples 
in the same way as the observations, and obtain values on 
900 grids. The best-fit parameter cr can be obtained from the 
grid which has the minimum value. As an illustration, in 
the lower panels of Fig. [I0]we compare the alignment signals 
obtained from the data (asterisks) with those obtained from 
the Monte Carlo simulations including the best-fit Gaussian 
distribution of misalignment angles (open symbols) . Clearly, 
the observed results can be well reproduced by such a model. 
Table 3 lists the best-fit values of a thus obtained, for a 
variety of cases that are shown in Figs. [6] - [9] To better 
describe these statistical values, 68.3% confidence levels are 
also given following each best fit cr. Here the confidence levels 
are obtained from the grids with = Xmm + 1-0. Clearly, 
the amount of scatter in the misalignment angle depends 
strongly on the group mass and on the color of the central 
galaxy. For the entire sample as a whole we obtain o ~ 
23°, while groups with red centrals, on average, have smaller 
misalignment angles than those with blue centrals. Owing 
to the nearly spherical shapes of dark matter haloes with 
12.0 < log[M//i"^M0] < 13.0, the misalignment angle could 
not be meaningfully constrained in these cases. 

In a recent study, K07 found that if the minor axis of 
the central galaxy and the angular momentum vector of the 
dark matter halo is in perfect alignment, most of the obser- 
vational satellite-central alignment signals can be automati- 
cally reproduced. On the other hand, any processes that can 
introduce an average misalignment ~ 40° between the minor 
axes of central galaxy and host dark matter halo can also 
explain the observed alignment signals. Using high resolu- 
tion numerical simulations, Bailin & Steinmetz (2005) found 
that the angular momentum vectors tend to align with the 
minor axes of the dark matter halo with a mean misalign- 
ment of ~ 25° at small radii and of ~ 40° at the halo virial 
radius. Here to be able to compare with our findings, as an 
experiment, we input misalignment angles (with an average 
~ 40° and some scatters) between the minor axes of the dark 
matter haloes and the central galaxies. By randomly project 
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Table 2. Best fitting parameters for the deviation angle between 
the orientation of the central galaxy and the dark matter halo. 



Subsample 


a (degree) 


alP 


90 q+O.l 

ZO.O_Q 3 


blue cen. 


38.ll°i 


red cen.*" 


16.6l°:l 


blue sat. 


29.3+E!-f 


red sat.'^ 


19.91-1 


red cen. /blue sat. 


38.5ii;:i 


red cen. /red sat. 


16 1+"--^ 


blue cen. /blue sat. 


80.01^1*^4 


blue cen./red sat. 




red cen. (M23<^) 




red cen. (M34) 


'•°-o.i 


red cen. (M45) 


25.71°:? 


blue cen. (M23) 




blue cen. (M34) 


64.81°:^ 


blue cen. (M45) 


9+0.1 


all (M23) 




all (M34) 


no Q + 0.2 
3 


all (M45) 


30.2l°i 



" all means all central galaxies 

blue cen. and red cen. denote blue centrals and red centrals, 
respectively. 

° blue sat. and red sat. denote blue satellites and red satellites, 
respectively. 

M23 means that the halo masses are in the range 
12 < log[M/ H-'^Mq] < 13. Note that in M23, the align- 
ment signals for the perfect alignment arc already as weak as the 
observed ones, and so no misalignment angle is introduced. 



the 3D shapes of the dark matter lialoes and central galax- 
ies, we measure the distribution of the angles between the 
projected major axes of the dark matter halos and central 
galaxies, if modelled with an Gaussian distribution, which 
roughly corresponds to a ~ 22°. Thus our finding are in re- 
markably good agreement with K07 and Bailin & Steinmetz 
(2005). 



5 SUMMARY AND DISCUSSION 

Using the large galaxy group catalogues constructed from 

the SDSS Data Release 4 (DR4) by Yang et al. (2007), 
we have investigated the shapes of their host dark matter 
haloes, and the correlation between the orientations of the 
central galaxies and those of their host haloes. In particu- 
lar, we obtained the two-dimensional and three-dimensional 
axis ratios of galaxy groups by comparing the observed, pro- 
jected elliptieity distributions of satellite galaxies with those 
of Monte Carlo simulations, and we determined the proba- 
bility distributions for the angles between the major axis of 
central galaxies and the lines connecting the centrals with 
their satellites. The main results of this paper are summa^ 
rized as follows. 

(i) Under the assumption that the spatial distribution of 
satellite galaxies traces the shapes of the underlying dark 
matter haloes we find that the projected elliptieity distribu- 



tions are slightly better fit with the triaxial models of JS02 
than with simple, spherical NFW models. 

(ii) The shapes of dark matter haloes depend strongly on 
their mass, with more massive haloes being more elongated. 
Haloes with masses in the range 12 < log[M/ /i~^Mq] < 13 
are nearly spherical, while more massive haloes with 14 < 
log[M/ h-^Me] < 15 are more prolate. 

(iii) There is no significant difference between the shapes 
of haloes with red central galaxies and those with blue cen- 
tral galaxies. 

(iv) Satellites are preferentially distributed along the ma- 
jor axes of their central galaxies. The strength of this 
alignment depends strongly on halo mass and on the col- 
ors of both central and satellite galaxies. The alignment 
is strongest between red centrals and red satellites, while 
blue centrals show almost no alignment at all. More massive 
groups show a stronger alignment than less massive groups. 

(v) The observed alignment can be reproduced if the pro- 
jected orientation of central galaxies is aligned with that of 
the projected mass distribution of its halo. However, the 
alignment is not perfect. The data can be reproduced under 
the assumption that the misalignment angle follows a Gaus- 
sian distribution around zero, and with a standard deviation 
of ~ 23 degrees. Because of (ii) and (iii) and (iv), groups 
with blue centrals have, on average, a larger misalignment 
angle than those with red centrals. 



Our findings regarding the shapes of the galaxy groups 
are in qualitative agreement with those of P06, and Plio- 
nis et al. (2004). However, there are quantitative difference, 
especially for small haloes. We find that small haloes are al- 
most spherical, and the axis ratios (major-to-minor) in our 
measurements are much smaller than their measurements. 
Note that we have used Monte Carlo simulations to infer the 
shapes of dark matter haloes, which avoids the impact of se- 
lection effects, while in Plionis et al.(2004) the axis ratios are 
directly measured from the distribution of member galax;ies. 
Although they found that poor groups are more elongated 
than rich ones, they cautioned that their results may be sig- 
nificantly affected by discreetness effects (see also Plionis et 
al. 2006). Numerical simulations (e.g, Kasun & Evrard 2005; 
AUgood et al. 2006) have shown that the small haloes are 
more spherical than massive haloes, in good agreement with 
our findings here. 

The alignment signals presented here are in good agree- 
ment with those obtained by Y06 using a similar, but smaller 
group catalogue constructed from the SDSS DR2, and with 
other studies (Brainerd 2005; Agustsson & Brainerd 2006a, 
2007; Donoso, O'Mill & Lambas 2006; Faltenbacher et 
al. 2007a; Azzaro et al. 2007) . Contrary to the previous stud- 
ies by AB06 and K07, our analysis of the implications for the 
correlation between the orientation of the central galaxies 
and that of its halo are based purely on the data presented 
here, and therefore does not depend on any galaxy formation 
models or numerical simulations. Nevertheless, tests show 
that their models of the perfect alignment between the mi- 
nor axis of the central galaxy and angular momentum vector 
of dark matter halo are in good agreement with our direct 
measurement of the projected misalignment angles. 
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